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ABSTRACT 

Three  different  strut  shapes  have  been  investigated  to 
determine  the  effect  of  shape  on  the  heave  motion  of  ships  at 
resonance.  A  time-domain  analysis  has  been  developed  to  include 
the  nonlinear  effect  of  large  motion  amplitude.  A  two-dimensional 
model,  in  which  the  ship  forward  speed  is  assumed  to  be  zero,  has 
been  investigated  to  simplify  the  problem.  The  results  show  that 
a  proper  choice  of  strut  shape  can  reduce  ship  heave  motion  sub¬ 
stantially  at  reonance. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  General  Hydromechanics  Research  Program, 
administered  by  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
(DTNSRDC)  Ship  Performance  Department,  and  was  authorized  by  the  Naval  Sea  Systems 
Command  (NAVSEA) ,  Hull  Research  and  Technology  Office.  Funding  was  provided  under 
Program  Element  61153N,  Task  Area  SR  0230101,  and  Work  Unit  1562-202. 

INTRODUCTION 

The  concept  of  a  small-waterplane-area  twin-hull  (SWATH)  ship  has  been  intro¬ 
duced  by  the  Navy  to  improve  the  seakeeping  characteristics  of  naval  ships  in  waves. 

A  SWATH  ship  consists  of  two  submerged  hulls,  and  struts  that  connect  the  wet  deck 

1* 

and  these  lower  hulls.  The  struts  are  usually  vertical-wall  sided. 

The  purpose  of  rhe  present  study  is  to  investigate  the  effect  of  strut  shape 
on  the  heave  motion  of  SWATH  ships  near  the  resonance  region.  In  particular,  we 
wish  to  determine  if  a  proper  selection  of  strut  configuration  will  decrease  the 
peak  value  of  heave.  An  existing  analytical  method  frequently  used  to  predict  the 
motions  of  SWATH  ships  is  based  on  a  frequency-domain  approach.  Since  the  motion 
amplitude  is  large  at  resonance,  the  frequency-domain  analysis  may  be  unable  to 
treat  the  nonlinear  effect  of  this  large-amplitude  motion.  In  this  study,  a  time- 
domain  analysis  has  been  developed  to  include  the  nonlinear  effect  of  large-motion 
amplitude.  To  simplify  the  problem,  a  two-dimensional  ship  model  in  which  the  ship 
forward  speed  is  assumed  to  be  zero,  has  been  investigated.  Extending  the  investi¬ 
gation  to  a  three-dimensional  SWATH  ship  is  beyond  the  scope  of  the  present  study 
and  should  be  considered  in  the  future. 


♦References  are  listed  on  page  21. 
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The  procedure  of  the  present  analysis  starts  with  computation  of  the  hydro- 
dynamic  and  hydrostatic  forces  for  three  different  ship  drafts  and  for  a  given  wave 
frequency.  The  equations  of  motion  are  integrated  numerically  at  each  time  step. 
Also,  at  each  time  step,  the  instantaneous  ship  draft  is  computed,  and  the  hydro- 
dynamic  and  hydrostatic  forces  are  interpolated  for  that  draft.  In  this  analysis 
it  is  assumed  that  the  wave  frequency  itself  is  constant  but  that  the  wave  amplitude 
increases  from  a  small  value  to  a  larger  steady-state  value.  This  study,  thus, 
attempts  to  analytically  duplicate  a  tank  test. 

Three  different  strut  shapes  have  been  investigated,  and  the  results  are  com¬ 
pared  with  those  obtained  using  the  frequency-domain  method.  The  computed  results 
show  that  a  proper  choice  of  strut  shape  can  reduce  the  heave  motion  results  sub¬ 
stantially  at  resonance. 

EQUATIONS  OF  MOTION 

We  define  an  oxy  coordinate  system  with  origin  0  at  the  free  surface.  The 
x  axis  is  in  the  plane  of  the  free  surface  and  the  y  axis  directed  vertically  upward 
along  the  SWATH  section  centerline  (see  Figure  1). 


Figure  1  -  Coordinate  System  Used  in  SWATH  Section  Study 

Let  us  assume  that  a  two-dimensonal  SWATH  section  is  floating  at  the  origin, 
and  that  a  train  of  waves  is  approaching  this  section  from  the  positive  x  axis. 
During  the  transient  period,  the  amplitude  of  the  wave  is  initially  very  small,  and 
it  increases  to  that  of  the  regular  wave  in  some  finite  time,  say  about  15  s. 
Furthermore,  it  is  assumed  that  the  wave  period  is  constant  during  the  transient 
period.  The  incoming  wave  sets  the  SWATH  section  into  a  harmonic  motion  and  the 
equation  for  the  heave  motion  can  be  written 
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(nri-A(t) )  y"  +  B(t)y'  +  C(t)y  +  D(t)y' 


=  F^(t)  cos  cot  +  F2(t)  sin  cot  (1) 

where  m  is  the  mass  of  the  section,  A  the  added  mass,  B  the  damping  force,  and 
and  F^  are  exciting  forces.  In  addition,  y  is  the  amplitude  of  heave  motion  and 
co  is  the  frequency  of  the  incoming  wave.  The  prime  and  double  prime  indicate  first 
and  second  derivatives  of  y  with  respect  to  time,  respectively.  The  time  dependency 

of  A,  B,  C,  D,  F. ,  and  F2  is  due  to  the  change  of  ship  draft  at  a  given  co  as  the 

wave  passes. 

The  following  computational  procedure  is  employed  for  obtaining  the  hydrostatic 
and  hydrodynamic  forces.  The  sectional  contour  is  divided  into  segments,  and  a  con¬ 
stant  source  of  unknown  strength  is  computed  using  the  body  boundary  condition 

which  states  that  the  normal  velocity  of  the  body  is  the  same  as  that  of  the  fluid 

2,3 

at  the  boundary.  This  method  is  sometimes  called  the  Frank-close-fit  method.  ’ 

If  we  assume  that  A,  B,  C,  F^,  and  F2  are  frequency  dependent  but  independent 

of  time,  and  that  D  is  zero,  then  the  solution  of  Equation  (1)  is  simply  in  the  form 

y  =  Re(y  e  (2) 


where  y  is  given  by 


y 


c 


+  iF 


2 


(m+A)  -  icuB 


(3) 


Equation  (3)  can  be  easily  obtained  by  substitution  of  Equation  (2)  into 
Equation  (1). 


NUMERICAL  INTEGRATION  OF  EQUATION  (1) 

When  the  coefficients  in  Equation  (1)  are  dependent  on  time,  it  is  impossible 
to  find  a  solution  such  as  Equation  (3).  The  solution  of  Equation  (1)  is  possible, 
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however,  through  a  step-by-step  numerical  integration.  The  Adams-Moulton  method  can 
be  used  to  integrate  Equation  (1)  numerically.  (The  details  of  this  method  are  given 
in  the  Appendix.)  In  order  to  utilize  the  Adams-Moulton  method,  we  rewrite  Equation 
(1)  as 


=  ay  +  by '  + 


cv'2  +  f. 


rne  Mf*  4-  f 


(A) 


where 


C 

m  +  A 


b 


B 

m  +  A 


D 

m  +  A 


f 


1 


m  +  A 


and 


m  +  A 


If  we  let 


dt  “ 


y' 


=  V, 


(5) 


then 
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With  v  =  z  these  equations  can  be  expressed  as  a  system  of  linear  equations 


and 


y'  =  z 

2 

z'  =  ay  +  bz  +  cz  +  cos  tot  +  f ^  sin  cot 


(7) 


With  substitution  of  x  for  t  we  can  now  apply  the  Adams-Moulton  method  explained 
in  the  Appendix  to  solve  Equation  (7). 

Before  we  start  the  numerical  integration  of  Equation  (7),  we  make  the  following 
assumptions : 

1.  The  displacement  and  KG  of  each  SWATH  section  are  constant. 

2.  The  draft  varies  as  the  wave  passes. 

3.  The  waterplane  and  beam  areas  at  the  free  surface  are  a  function  of  the 
draf  t . 

4.  The  added  mass,  and  damping  and  exciting  forces  are  a  function  of  the 
draft  at  a  given  frequency  co. 

In  the  present  method,  the  coefficients  A,  B,  C,  F^,  and  in  Equation  (]) 
are  computed  for  three  different  drafts  at  a  given  frequency  using  the  frequency- 
domain  ship  motion  program;  one  of  these  is  the  calm-water  draft.  At  the  initial 
stage,  the  values  of  these  coefficients  are  for  the  calm-water  draft.  When  the 
incoming  wave  approaches  the  body,  an  instantaneous  draft  is  computed,  and  those 
coefficients  are  interpolated  for  that  draft. 

NUMERICAL  RESULTS 

Three  two-dimensional  SWATH-ship- like  sections  (Figure  2)  have  been  selected 
as  a  numerical  tesl  in  the  present  study.  SWATH  Section  A  is  almost  the  same  as 
the  midship  section  of  one  hull  of  SWATH  6A.^  The  strut  of  this  section  is 
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At  the  start  of  computation  (or  n  =  o)  the  values  of  the  backward  difference  are 
not  available.  We  apply  the  Taylor  series  method  to  calculate  these  tabular  values. 
The  Taylor  series  for  y  is  given  by 


y  = 


y  +  y ' (x-x  )  + 
o  o  o 


2! 


<x-xo)‘ 


3! 


<x-Xo} 


(A5) 


where  y^  is  the  value  of  Equation  (Al)  at  x  =  x^,  y^'  the  value  of  y'1  at  x  =  x  , 
and  so  on.  y^'  can  be  obtained  after  differentiating  Equation  (Al)  and  substituting 
x  =  x^,  y  =  yQ,  and  y'  =  y^.  In  order  to  compute  y  or  yn  ^  for  n  =  0,  we  compute 
y'o,  y'Q' ,  y^'  ' ,  ...  with  initial  values,  and  then  calculate  y  by  means  of  Equa¬ 
tion  (A5)  with  an  assumed  value  for  (x  ,  -  x  ) .  With  a  known  value  of  y  ,  we 

-I  o  -I 

calculate  the  derivatives  y^,  y^«  y '  j '  ,  ...  ,  and  y_^  can  be  calculated  with 
substitution  of  (x  ^  ~  x_^)  for  (x  ^  -  XQ)  •  Once  y  y'^,  ^-2’  '  '  ‘  are  known,  t*ie 
backward  differences  of  y'  are  easily  computed. 

When  a  system  of  first-order  differential  equations  is  given  as 

y’  =  f (x,y , z) 
z'  =  g(x,y,z) 

with  initial  values  of  x  =  x  ,  y  =  y  ,  and  z  =  z  ,  the  solution  can  be  obtained  by 
the  Adams-Moulton  method.  In  this  case  we  need  two  tables  of  backward  differences, 
one  each  for  y  and  z.  Equation  (A5)  is  used  to  calculate  the  values  at  the  start 
of  the  computation  for  y  and  z.  Once  the  tables  are  ready,  y^+^  and  zn+^  are 
computed  with  an  equation  having  the  form  of  Equation  (A2) ,  and  we  call  these  the 
first  results.  A  new  table  of  backward  differences  is  computed  with  the  first 
results.  That  is,  the  second  results  for  y  ,,  and  z  ,,  are  computed  with  Equation 
(A2)  and  this  yields  a  new  table  of  backward  differences.  These  values  are  compared 
with  the  first  results.  If  the  convergence  is  not  within  the  desired  limit,  the 
procedure  is  repeated  using  the  second  results  to  compute  third  results,  and  so  on, 
until  the  convergence  is  satisfactory. 
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Vy'  V2y'  V3y'  V4y’ 


n-3 


r 

n-3 


n-3 


v2c2 


7\-l 


Vy- 


n-2 


,3Ci 


n-2 


n-2 


n-2 


y2Ci 


V4y' 

n 


yCi 


y3*; 


n- 1 


n-1 


n- 1 


The  numerical  procedure  of  the  Adams  method*  is 

1.  Use  Equation  (A2)  to  compute  y^+^ 

2.  From  the  given  differential  equation  Equation  (Al)  ,  calculate  y'n+^ 

3.  Complete  the  difference  table  for  the  (n+l)jdt  line 

4.  Using  the  (n+l)^h  line  in  place  of  the  nth  line,  repeat  the  three  steps. 
Equation  (A2)  computes  Yn+^  without  using  y^+^.  Moulton’s  method,  however, 

computes  Yn+^  using  y^+^  more  accurately.  The  solution  of  Equation  (Al)  by 
Moulton's  method  is  given  by 


*n+i  =  *n  +  hfy;+r<1/2>vy;+r(1/24>vX+i 


(19/270)V4y^+1-(27/1440)V5y^+1+  .. .] 


(A4 ) 


♦Details  of  this  method  are  given  in  W.E.  Milne,  "Numerical  Solution  of  Differ¬ 
ential  Equations,"  Dover  Publications,  Inc.,  New  York  (1970). 
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APPENDIX 

ADAMS-MOULTON  METHOD 


When  a  differential  equation  is  given  as 


f(x,y) 


(AD 


with  initial  values  x  =  xq  and  y  =  yQ,  the  approximate  solution  can  be  expressed 
with  the  backward  difference  as 


y  =y  +  h[y'+(l/2)Vy,+(5/12)V2y'+(3/8)V3y'+(251/720)V4y,+...]  (A2) 

n+i  nn  n  n  n  n 

where  h  is  the  step  increase  of  x,  and  the  backward  difference  operator  V  is  defined 
as 


V y'  =  y '  -  y  , 
n  yn  'n-1 


V2y'  =  Vy'  -  Vy'  .  =  y'  -  2y'  .  +  yr  „ 
7n  ^n  Jn~l  n  7n-l  n-2 


V 3y '  =  V2y'  -  V2y'  .  =  y '  -  3y '  .  +  3y '  -  y ' 

Jn  n  yn-l  n  yn-l  n-2  Jn-3 


(A3) 


and  so  on.  The  computational  scheme  of  the  backward  difference  is  given  in  tabular 
form  as 


As  shown  in  Figure  11,  the  heave  characteristic  for  Section  C  is  opposite 
that  of  Section  B  in  that  positive  amplitude  is  generally  smaller  than  the  negative 
one.  In  addition,  the  Section  C  hydrostatic  and  exciting  force  coefficients  also 
exhibit  a  different  trend. 

Tible  2  shows  that  if  we  introduce  viscous  damping  (C^  =  0.715  instead  of  0), 
a  further  reduction  of  heave  motion  is  obtained. 

SUMMARY  AND  CONCLUSIONS 

Time-domain  analysis  is  applied  to  investigate  the  effect  of  SWATH  ship  strut 
shape  on  heave  motion  at  heave  resonance.  The  computed  results  show  a  strong  depen¬ 
dence  of  the  heave  motion  on  strut  shape.  A  substantial  reduction  in  motion  is 
achievable  by  proper  strut  selection  for  SWATH  ships.  The  following  conclusions 
may  be  drawn  from  the  present  study: 

1.  Heave  motion  at  resonance  can  be  reduced  substantially  by  a  proper  choice 
of  strut  shape.  Furthermore,  there  is  an  additional  reduction  of  the  ship  motion 
amplitude  when  viscous  damping  forces  are  incorporated. 

2.  The  present  study  was  carried  out  only  for  two-dimensional  sections.  It 
is  of  interest  to  conduct  a  similar  investigation  for  existing  SWATH  ship  models, 
such  as  SWATH  6A,  6B,  6C,  and  6D.  We  can  then  compare  the  results  with  experimental 
data  obtained  in  regular  waves. 
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The  motion  results  for  variable  draft  for  SWATH  Sections  B  and  C  at  resonance 
are  substantially  reduced  compared  to  Section  A.  Furthermore,  for  the  present 
method,  the  location  of  resonance  is  slightly  shifted  towards  the  lower  frequencies 
(see  Figures  3(b)  and  (c)).  As  shown  in  Figures  4  and  5,  the  added  mass,  damping, 
hydrostatic  and  exciting-force  coefficients  for  Sections  B  and  C  vary  significantly 
with  draft  change. 

Figures  6  to  11  show  heave  motion  time  histories  for  the  three  sections  at 
resonance  with  =  0.  For  each  SWATH  section,  the  motions  for  constant  dratt  uu.l 
variable  draft  are  computed.  For  the  case  of  constant  draft,  the  values  of  added 
mass,  damping,  exciting,  and  hydrostatic  forces  are  taken  when  the  draft  h  is  8.1  m; 
this  is  true  throughout  the  total  time  of  400  s.  As  mentioned  earlier,  the  result 
of  the  present  method  when  assuming  constant  draft  should  be  the  same  as  those  of 
the  frequency-domain  approach.  Figures  6  to  8  show  the  resul  s  for  constant  drat i 
The  heave  amplitudes  have  not  reached  a  steady-state  value  by  t  =  400  s.  When  t  i  . 
taken  as  1000  s  (which  is  not  shown  in  these  plots),  heave  becomes  steady  state,  a.,  t 
the  amplitudes  are  then  the  same  as  the  frequency-domain  results.  Most  runs  were  n.  i 
made  for  1000  s,  however,  because  of  excessive  computer  time  and  cost. 

Figures  9  to  11  show  the  results  for  variable  draft.  Since  SWATH  Section  A  has 
a  vertical-wall  sided  strut,  there  is  no  change  in  the  heave  time  history  compared 
to  the  constant-draft  case.  For  SWATH  Sections  B  and  C,  however,  the  amplitudes 
at  resonance  become  much  smaller  than  those  for  constant  draft.  The  rapid  reduction 
of  the  resonance  value  of  heave  motion  for  SWATH  Sections  B  and  C  is  caused  by  the 
fact  that  the  hydrostatic  and  hydrodynamic  force  coefficients  (Figures  4  and  3) 
change  greatly  as  a  function  of  draft.  These  same  coefficients  for  SWATH  Section  A 
are  less  dependent  on  the  draft. 

It  is  worth  noting  that  the  heave  amplitudes  for  variable  draft  for  SWATH 
Sections  B  and  C  are  not  symmetric  with  respect  to  the  undisturbed  free  surface. 

For  SWATH  Section  B,  the  positive  amplitude  (heave  up)  tends  to  be  larger  than  the 
negative  one  as  the  motion  continues  (Figure  10).  The  hydrostatic  and  hydrodynami. 
force  coefficients  for  Section  B,  conversely,  become  smaller  for  upward  heave 
(Figure  4).  Thus,  even  though  the  motion  is  a  function  of  the  force  coefficients 
according  to  Equation  (3),  it  appears  that  small  coefficients  result  in  large  heave, 
and  vice  versa. 
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TABLE  1  -  PRINCIPAL  DIMENSIONS  OF  SWATH  SECTIONS 


Principal 

Dimensions 

(m) 

SWATH  Section 

A 

B 

C 

a 

2.3 

2.3 

2.3 

h 

8.1 

8.1 

8.1 

f 

5.8 

5.8 

5.8 

i 

2.02 

2.16 

1.90 

bl 

1.1 

1.41 

0.45 

b2 

- 

0.4 

1.3 

KG 

3.01 

3.01 

3.01 

S* 

24.52 

24.52 

24.52 

2 

*Sectional  area  (m  ) . 

TABLE  2  -  NONDIMENSIONAL  AMPLITUDE  (£/A')  AT  HEAVE  RESONANCE 


SWATH  Section 

A 

B 

C 

Resonant  Frequency 
oj  (rad/s) 

0.780 

0.855 

0.705 

o 

c  a 

0.0 

/ 17.56 

13.13 

■U  AJ 
(0  M-l 

e  to 
o  u 
o  a 

,  0.715 

< 

\  / 

3.53 

3.31 

<u  o 
«h  a 

X 3 

0.0 

W' 

17.52 

3.24 

2.23 

(t!  -U 
■H  U-l 

m  m 
nl 

>  n 

,  0.715 

3.53 

2.79 

2.11 

Frequency-Domain 

Method 

18.84 

13.64 

32.78 

vertical-wall-sided.  SWATH  Section  B  has  a  flare-type  strut  with  minimum  width 
at  the  intersection  with  the  lower  hull.  SWATH  Section  C  has  a  maximum  width  at 
the  intersection  with  the  lower  hull  and  a  minimum  width  at  the  free  surface.  All 
three  sections  have  the  same  radius  a  for  the  submerged  circular  hull.  The  principal 
dimensions  of  these  sections  are  indicated  in  Figure  2  and  tabulated  in  Table  1. 

The  numerical  computation  has  been  performed  for  two  cases:  (a)  constant  drat  t 
(the  results  of  this  case  are  the  same  as  for  those  of  the  f requncy-domain  computa¬ 
tion),  and  (b)  the  draft  varies  as  the  wave  passes.  For  each  case  the  computation 
was  done  with  and  without  inclusion  of  the  viscous  damping  force  D.  The  viscous 
damping  force  is  computed  as 

D  =  (pCDsv2)/2 

where  s  is  the  waterplane  area  and  v  is  the  instantaneous  vertical  velocity  of  the 
SWATH  section.  Whenever  this  damping  force  was  included  in  the  computation,  it  was 
assumed  that  the  coefficient  CQ  is  0.715. 

Figure  3  shows  the  heave  motion  results  for  the  three  sections  selected.  These 
results  are  compared  with  responses  predicted  by  the  frequency-domain  motion  program. 
Since  SWATH  Section  A  has  a  vertical  strut,  there  is  no  change  in  the  response 
function. 

In  the  plots  for  SWATH  Sections  B  and  C  the  peak  values  obtained  by  time-domain 
analysis  (present  method)  are  substantially  smaller  than  the  value  for  Section  A. 

This  is  true  because  the  struts  of  Sections  B  and  C  provide  greater  damping  than  the 
wall-sided  strut  does.  Table  2  shows  the  results  of  heave  motion  at  resonance.  The 
results  obtained  for  constant  draft  with  =  0  should  be  the  same  as  those  of  the 
frequency-domain  computation.  However,  there  is  a  substantial  difference  in  the 
results  for  SWATH  Section  C.  This  difference  is  caused  by  the  fact  that  in  the 
present  method  the  total  integration  time  of  the  differential  equation  is  of  too 
short  a  duration.  Throughout  the  computations  this  time  was  taken  as  400  s.  If  a 
total  integration  time  of  1000  s  is  used,  the  results  obtained  for  SWATH  Section  C 
are  the  same  for  both  the  time-domain  and  frequency-domain  computations. 
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